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Abstract. I have searched the Westerbork Northern Sky 
Survey (WENSS) source list for detections of known ra- 
dio pulsars. A source with a flux density greater than five 
times the local noise level is found near the positions of 25 
pulsars. The probability that one out of these 25 sources 
is a chance coincidence is about 10%. I have looked at the 
WENSS maps of the non-detected pulsars. A flux den- 
sity between three and five times the local noise level is 
found near the positions of 14 of these non-detected pul- 
sars. There is a 50 percent probability that (at least) one of 
these marginal detections is just a noise fluctuation. Four- 
teen radio pulsars, which according to earlier flux measure- 
ments have flux densities above three times the WENSS 
noise level, are not detected. Of the 39 pulsars detected in 
the WENSS 19 are also detected in the NRAO VLA Sky 
Survey (NVSS). By combining the WENSS and NVSS flux 
densities for these 19 pulsars spectral indices are obtained 
that differ by up to 50% from the long term averaged val- 
ues reported in the literature. This affects the reliability 
of pulsar candidates that are selected on the basis of their 
WENSS-NVSS spectral index. 
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The pulsar population found in this way, may supple- 
ment the presently known population, since this method 
has totally different selection effects. These effects can be 
investigated by studying the spectral indices and polari- 
sation degrees of known pulsars in these continuum ob- 
servations. The spectral indices are strongly influenced by 
scintillation. This may cause the flux density of a pulsar to 
vary by more than 100 percent on time scales of minutes 
(diffractive scintillation) to days (refractive scintillation). 

In this paper I search for detections of pulsars in the 
Westerbork Northern Sky Survey (WENSS), a survey per- 
formed at 325 MHz . I compare my res ults w ith those from 
Kaplan et al. (|l998|) and Han & Tian (|l999| ), who did sim- 
ilar analyses with data from the 1400 MHz NRAO VLA 
Sky Survey (NVSS). Sect. | describes the WENSS and 
Sect. |^ describes how the pulsar catalog was correlated 
with the WENSS source catalog. In Sect. |] the positions 
of the non-detected pulsars are searched for flux densities 
above three times the local noise level. Sect. [| discusses the 
remaining non-correlations. Sect. || combines the WENSS 
results with the NVSS correlation studies to determine 
the spectral indices and compares these with values re- 
ported in other literature. Finally, in Sect. ^ the role of 
scintillation is discussed. 



1. Introduction 

Pulsars can have a steep radio spectrum at frequencies 
around 1 GHz and can be highly polarised (e.g. Manch- 
ester & Taylor 1977) . Synthesis maps can therefore be 
used to select pulsar candidates on the basis of their 
steep spectrum and/or high degree of polarisation. These 
sources are later observed with a high time resolution in- 
strument to search for pulsations. At least two pulsars, 
PSR B1937+21 (packer et al., 1982[ ) and PSR J0218+4232 
(Navarro et al., 1995), have been found in this way, while 



they were missed in regular pulsar surveys because their 
pulses were smeared out in the detection process due to 
their small period and their high dispersion measure. 
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2. WENSS 

The Westerbork Synthesis Radio Telescope (WSRT) is 
an east-west array, consisting of fourteen 25 meter dishes 



(Baars & Hooghoudt, 1974). A twelve hour observation is 



needed to get spatial resolution in all directions. Ten of 
the dishes have fixed positions and are spaced 144 meters 
apart. The remaining four are moveable, although their 
mutual distances are usually kept constant. By observing 
an object several times with different distances between 
the fixed and moveable subarray, different baselines are 
obtained and the uw-coverage improves, yielding an im- 
proved synthesised antenna pattern. 

The Westerbork Northern Sky Survey (WENSS) is a 
325 MHz continuum survey of the sky above declination 
+30° (Rengelink et al. 1997, de Bruyn et al. in prepara- 
tion). This area was surveyed using a mosaicing technique. 
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nation > +75°) was surveyed with a larger bandwidth 
and the detection limit for this area was about 10 - 12 
mJy. The WENSS source catalog contains 229420 sources 
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Fig. 1. Pulsar positions and WENSS source densities. 
Filled circles mark the positions of the pulsars which are 
correlated with a WENSS source. Open circles indicate 
the pulsars, with a marginal detection in the WENSS. 
Confused pulsars are marked with a star. Pulsars with no 
WENSS counterpart are indicated with a plus. The grey- 
scale indicates the WENSS source density (white is low 
density). The polar cap region has been surveyed with a 
larger bandwidth and has a higher source density. A small 
region near Cygnus A (RA 20h, Dec 40°) is blanked in the 
WENSS. The solid line indicates the Galactic plane. 

Each field was observed 18 times for 20 seconds spread 
over 12 hours. Each mosaic was observed on six days with 
a different spacing between the fixed and moveable sub- 
array to get a uniform spatial distribution. These obser- 
vations were spread over periods of weeks to years. The 
resulting flux densities are averaged over all these obser- 
vations. 

The WENSS beam size is 54" x 54" cosec S (FWHM). 
The final maps have pixel sizes of 21.09". When a pixel 
was found with a flux density above five times the local 
noise level, a two-dimensional Gaussian was fitted to its 
surroundings. The coordinates of the centroid of the fit, 
the maximum (peak) flux density and the flux density in- 
tegrated over the fit were added to the source list. Extra 
flags in the catalog mark multiple and extended sources. 
For a point source the peak flux density equals the flux 
density integrated over the beam. Only the peak flux den- 
sities are used in this analysis, since pulsars are intrinsi- 
cally point sources. Bright sources have positional errors 
of 1.5", weaker sources at lower declinations have errors 
up to about 10". On average the uncertainty is 5". 

The total bandwidth was 5 MHz. For most regions of 
the sky the detection limit was between 15 and 25 mJy 
(five times the local noise level). The polar cap (decli- 



( |dc Bruyn et al., 1998[ ). A total of 18186 of these are lo- 
cated in the polar cap area. Fig. [j] shows the density of 
sources in the WENSS area. 



3. Detected pulsars 

3.1. Positional correlation 

The radio pulsar catalog (Taylor et al., 1993, Taylor et 
al., 1995) contains 84 entries in the part of the sky that 
was covered by the WENSS area. The typical positional 
uncertainty for a pulsar with a flux density greater than 10 
mJy is 0.1" or less. In most cases the uncertainty in the 
pulsar position is negligible compared to the positional 
uncertainty of WENSS sources. 

The pulsar proper motions are neglected, since for each 
pulsar the change in position between the epoch of discov- 
ery and the epoch of the WENSS is less than 0.5" in right 
ascension and less than 0.8" in declination for all these 
pulsars. This is much smaller than the uncertainties in 
the WENSS positions. 

Seven pulsars in the WENSS area have large posi- 
tional errors of about 4'. The probability that a WENSS 
source is located by chance within a circle of three times 
this positional error is more than 90 percent, if it as- 
sumed that the WENSS sources are uniformly distributed 
over the sky. Therefore, I have excluded these 7 pul- 
sars (PSRs J0417+35, B1639+36B, J1758+30, J1900+30, 
J1931+30, J2002+30 and J2304+60) from further analy- 
sis. PSR B2000+40 was also excluded, since it is located 
in the Cygnus A region, where no WENSS map could be 
made. 

I have taken the J2000 positions of the remaining 
76 pulsars and compared them with the positions of the 
sources in the WENSS catalog. Twenty-five pulsars have a 
WENSS source located within three times their combined 
positional uncertainty a, with 




(1) 



where Aa and AS are the positional differences in right 
ascension and declination, respectively, and 



'a. WENSS + °a,PSR (^) 

and similarly for as- Table [j] lists the J2000 positions of the 
correlated pulsars and WENSS sources and their offset, 
both in arcseconds and in a. 

In Fig. [I] it can be seen that the pulsars are located in 
areas with different values for the WENSS source density. 
The probability of a change coincidence can be approxi- 
mated by a probability calculation that assumes a uniform 
distribution of the WENSS sources. In that case, the prob- 
ability that an individual correlation is just by chance is 
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Fig. 2. Histogram of differences between the position of 
a pulsar and its nearest WENSS source in units of the 
combined positional uncertainty a. The dotted line markes 
the 3er limit. The dashed lines is the expected distribution 
for 25 correlations, i.e. the number of found pairs with a 
positional difference below 3c. 



0.0012. The binomial probability, that one out of 76 trials 
gives a chance correlation is 0.083. The probability that 
two correlations occur by chance is 0.004. 

Fig. ^displays the distribution of positional differences 
between a pulsar and its nearest WENSS source in units of 
their combined positional uncertainty a. There is a clear 
gap between the correlated pairs (difference less than 3u) 
and the non-correlated ones. The distribution of the posi- 
tional difference (A) for the related pairs is 



P(A) = A 



-A-72 



(3) 



which follows after a conversion to polar coordinates. 
This distribution is also plotted in Fig. |^. A Kolmogorov- 
Smirnov test assigns a 15 procent probability that our 
sample is drawn from the distribution (j^) . 

There are two objects with positional differences be- 
tween 3 and 10 a. These are in confused regions and will 
be discussed in Sect. [|. 

3.2. Flux densities 

Table || lists the observed flux densities of the correlated 
WENSS sources and their uncertainties. The flux densi- 
ties can be compared with known pulsar flux densities, 
but these are usually measured at frequencies of 400 MHz 
and higher. By assuming a power law with constant spec- 
tral index, these flux density data can be extrapolated. I 
have used flux density data from Lorimer et al. ( 1995| ), 
hereafter LYLG. These flux densities are averaged over 
many observations spread over years and their uncertain- 
ties include variations due to scintillation. LYLG provide 



Table 2. List of estimated pulsar flux densities (PSR 
S'325) and WENSS source peak flux densities (WENSS 
S'325) for the 25 correlated objects. Navarro et al. 
(1995) observed that the continuum flux density for 
PSR J0218+4232 varies significantly between 100 and 200 
mjy. 
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data for 24 of the 25 pulsars in Table [l]. Flux densities 
for PSR J0218+4232 are taken from its discovery paper 
( [Navarro et al., 1995| ). 



Pulsar flux densities usually obey a power law with a 
negative exponent (5 oc v a with a < 0) in the frequency 
range from 325 to 1400 MHz. I have fitted the logarithms 
of the flux densities with a straight line. These lines are 
plotted in Fig. ||. From this fit a flux density at 325 MHz 
is estimated. This estimate is plotted against the flux den- 
sity of the WENSS counterpart in Fig. ^. It is known that 
some pulsars have a low frequency turnover, usually lo- 
cated around 100 MHz (|Malofeev et al., 1994]). However, 



some pulsars exhibit a turnover at a higher frequency : 
PSR B0329+54 around 300 MHz ( |Lync fc Rickett, 1968|) 
and PSR B2021+51 around 400 MHz. The spectrum of 
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Fig. 3. Flux densities against frequency for 24 pulsars with a WENSS counterpart (PSR J0218+4232 is excluded, 
see text). Filled circles mark flux densities measured by Lorimer et al. (1995). The dotted line is a fitted power law 
with a constant spectral index. Flux densities from the corresponding WENSS sources are indicated by a cross. Flux 
densities from the corresponding NVSS sources are marked with a plus. The errors in the WENSS and NVSS source 
flux densities are smaller than or of the order of the size of the symbol and do not include any uncertainty caused by 
scintillation. 
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Table 1. J2000 positions of pulsars and correlated WENSS sources. Offsets are the roots of the sum of the squared 
differences in right ascension and declination in " (absolute) and in a (relative). 
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PSR 1 32319+60 is flat between 200 and 600 MHz (|Mal- 
ofeev e jt al., 1994[) . PSRs B1946+35 and B2154+40 also 
have a turnover, but it is not clear whether this located 



between 325 and 400 MHz ( |Malofeev et al., 1994|) . In these 
cases the assumption of a constant power law is not cor- 
rect and this results in an overestimated flux density at 
325 MHz. 



Navarro et al. (1995) discovered that the flux density 
of PSR J0218+4232 has a non-pulsed component. They 
find that the continuum flux density at 325 MHz varies 
between 100 and 200 mJy. I used this flux density esti- 
mate for a comparison with the WENSS flux density. An 
estimate from the pulsed flux density was not derived. 

3.3. Source maps 

Fig. U shows maps of the WENSS sources which are cor- 
related with a pulsar. Most sources are consistent with a 
point source convolved with the beam shape. The shapes 



of PSR B0329+54 and PSR B0809+74 are affected by 
scintillation (see Sect. 0). The integrated and peak flux 
densities of PSR J0218+4232 differ by 20 percent. This 
suggests that the source is extended. However, a two- 
dimensional fit shows that its shape is not significantly 
different from the beam shape. 



4. Marginally detected pulsars 

The WENSS source finding routines have a flux density 
threshold at five times the local noise level (SS'noisc)- Since 
in this paper single trials are done to find a correlation at 
a very small number of positions, a 3S , no i S o source can be 
marked as a marginal detection. However, source fitting 
is less acurate. Therefore, a 5 by 5 pixel box (1.75') cen- 
tered around the pulsar position is searched for the pixel 
with the maximum flux density. This pixel was marked 
as a marginal detection, if its flux density was greater 
than SS'noisc. The positional uncertainty for these marginal 
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Fig. 5. Maps of pulsars that are detected in the WENSS. The plots are centered around the pulsar position, marked 
with a cross. The grey oval in the lower left corner indicates the FWHM beam size. Contour levels are at 5, 10, 15, 
20, 28, 40, 64, 128, 256, 512 and 1024 mJy (drawn lines) and at -5 and -10 mJy (dotted lines). 



sources is estimated using the same equation as for the fit- 
ted sources (Rengclink ct al., 1997). The positional error 
for these marginal sources is approximately 42". Since a 
1.75' box is searched, no marginal detection is missed. 

The probability of finding a 3S no i sc (or higher) pixel in 
a map with Gaussian noise is 2.7 • 10~ 3 . About 5 WENSS 
FWHM beams fit into a 5 x 5 pixel box. Fifty-one boxes 
were searched. Effectively, 255 trials have been performed. 



The binomial probability that one 3S no i sc pixel is found is 
0.35, the probability that two are found is 0.12 and that 
three are found is 0.03. Therefore, there is a 50 percent 
probability that (at least) one of our 14 marginal detection 
is just a noise fluctuation. 

The maps of the marginal detections are shown in Fig. 
[| The location of the pixel with maximum flux density 
is also indicated. The shapes of these sources is not as 
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Fig. 6. WENSS maps of the marginal detections which are correlated with a pulsar. The plots are centered around the 
pulsar position, marked with a cross. The dots indicate the maximum pixel in a 5 x 5 pixel box around the centrum. 
The grey oval in the lower left corner indicates the beam size. Contour levels at 5, 10, 15, 20, 28, 40, 64, 128, 256, 512 
and 1024 mJy (drawn lines) and at —5 and —10 mJy (dotted lines). 



point-like as the strong detections with S > 5S nolsc . The 
flux densities of the marginal detections are listed in Table 
together with an estimate of the pulsar flux density 
based on a simil ar e xtrapolation of the pulsar spectrum 
as done in Sect. |3~2|. The estimate for PSR J1518+4904 



is base d on its measured flux density at 370 MHz ( Bayer 
et al., 1997), since its spectrum as plotted by Kramer et 



al. ( 1999| ) shows evidence for a low frequency turnover. 

The ratios of the extrapolated pulsar flux density and 
the WENSS source flux density are displayed in Fig. ||. 
The spread is of the order of a factor 1.5, which is compa- 
rable with the spread for the detected sources that were 
discussed in the previous section. Five sources were de- 
tected, although their expected flux density was below 
three times the local noise level (see Table ||). The flux 
density at 325 MHz could not be estimated for three other 
pulsars, since no reliable flux density data at other fre- 
quencies were available. 

The contours of four pulsars are confused by nearby 
radio sources. These sources are discussed in the following 
and are shown in Fig. [7]. 

PSR B0655+64: the WENSS source and the pulsar 
position are 4.3cx apart. The estimated flux density of the 



pulsar at 325 MHz is 7 ± 2 mJy, but the WENSS source 
is 23 ± 5 mJy. Also, the NVSS (see Sect. 0) shows a radio 
source at the WENSS position and clearly away from the 
pulsar position. Its flux density at 1400 MHz is about 5.6 
mJy, while the pulsar flux density is expected to be 0.3 ± 
0.1 mJy. The pulsar has a proper motion, but it is small 
and directed towards negative declinations. 

PSR B1112+50: The WENSS source is bright (135 
mJy) and has an accurate position. The separation be- 
tween the fitted WENSS position and the known pulsar 
position is 4.6c. Extrapolation of the pulsar spectrum re- 
sults in an estimated flux density at 325 MHz of 16 ± 3 
mJy, much less than that of th e WE NSS source. Kaplan 
et al. ([L998D and Han & Tian (|1999| ) searched the NVSS 
for pulsar counterparts and also noted that the pulsar is 
confused by a strong NVSS source 12" away. 

PSR B1951+32: The coincident source in the WENSS 
catalog is marked as extended and the emission is dom- 
inated by the supernova remnant CTB80. The pulsar is 
associated with this remnant (Strom 1987, Kulkarni et al. 
1988). The WENSS peak flux density is 983 mJy, which is 
about a factor 70 stronger than the expected pulsar flux 
density. 
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Fig. 7. WENSS maps of the pulsars, which are confused by other radio sources. The plots are centered around the 
pulsar position, marked with a cross. The grey oval in the lower left corner indicates the beam size. Contour levels at 
5, 10, 15, 20, 28, 40, 64, 128, 256, 512 and 1024 mJy (drawn lines) and at -5 and -10 mJy (dotted lines). 




10 100 

PSR S 325 (mJy) 



1000 



Fig. 4. WENSS flux density as a function of the extrap- 
olated pulsar flux density (PSR S325) for the 25 corre- 
lated objects (filled symbols) and for the marginal ob- 
jects (clear symbols). PSR J0218+4232 is marked with a 
square. PSR B0329+54 and PSR B2021+51 show indica- 
tions for a low frequency turnover and are marked with a 
triangle. 



PSR B2306+55: It can be clearly seen in the map that 
this source has two components, of which the weaker one is 
probably the counterpart to the pulsar. This component is 
not listed in the WENSS source list. The estimated pulsar 
flux density is about 30 ± 3 mJy. The bright component of 
the WENSS source is 125 ± 9 mJy, separated 28 a from the 
pulsar position. The second component has a flux density 
of approximately 24 mJy. 



Table 3. List of estimated pulsar flux densities (PSR 
S325), WENSS maximum pixel flux densities (WENSS 
S'325), the local noise level in the WENSS map 

(»5noise) 

and the ratio R of the WENSS flux density and the noise 
level for the 14 marginal detections. PSR B0053+47 and 
PSR B0153+39 have only upper limits for their spectral 
index, so they have a lower limit for their 325 MHz flux 
density. The estimate for PSR J1518+4904 is based on its 
370 MHz flux density (see text). 'D' in the last column 
indicates that a source near the pulsar position was de- 
tected, although its expected flux density was below three 
times the local noise level. 



PSR 
name 


PSR 5*325 

estimate 
(mJy) 


WENSS 

5325 
(mJy) 


(mJy) 


R 






B0037+56 


10.3 


± 


0.7 


11 


± 


5 


3.6 


3.2 


D 




B0053+47 


> 4.2 


± 


0.6 


13 


± 


5 


3.5 


3.7 


D 


? 


B0105+68 


5.6 


± 


0.8 


16 


± 


5 


3.5 


4.6 


D 




B0153+39 


> 6.7 


± 


1.8 


12 


± 


5 


3.5 


3.4 


D 


? 


B0154+61 


8.4 


± 


1.1 


14 


± 


5 


4.0 


3.4 


D 




B0353+52 


19 


± 


3 


14 


± 


5 


3.9 


3.5 






B0458+46 


18.4 


± 


1.0 


18 


± 


6 


4.6 


4.0 






B0609+37 


21 


± 


3 


13 


± 


5 


3.8 


3.3 






J1518+4904 


~ 8 


± 


4 


13 


± 


4 


3.0 


4.2 


D 


? 


B1905+39 


40 


± 


5 


16 


± 


5 


3.8 


4.2 






B1953+50 


37 


± 


6 


29 


± 


9 


6.5 


4.4 






B2227+61 


28.6 


± 


1.7 


18 


± 


6 


4.7 


3.8 






B2323+63 


11 


± 


3 


21 


± 


8 


6.0 


3.5 


D 




B2324+60 


21.7 


± 


1.8 


27 


± 


12 


8.9 


3.0 


D 





In all these four cases I conclude that the source in the 
WE NSS catalog and the pulsar are unrelated. Galama et 
al. ( 1997 ) reached the same conclusion for PSR B0655+64. 
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5. Non-detected pulsars 

Table f| lists the pulsars that have no counterpart in the 
WENSS. In 14 cases the expected pulsar flux density is 
higher than three times the local noise level. Still, the pul- 
sar was not detected. In two cases no reliable pulsar flux 
density estimate at 325 MHz is available. In three other 
cases the estimate is based on 400 MHz observations. In 
case of PSRs B0841+80 and B1839+36A this was done, 
because there was no spectral information available. The 
spectrum of PSR J1012+5307 might also have a low fre- 
quency turnover (see its spectrum as plotted by Kramer 
et al. 1999) . Its flux density is known to vary by up to a 



factor four from its mean value of 30 mJy (Nicastro et al 
19950 

Five (and possibly eight) pulsars are detected, which 
were expected to be not detectable. The number of non- 
detected sources that were expected to have a flux density 
greater than the detection limit, should be roughly the 
same as the number of unexpected detections. The dif- 
ference may be due to Poisson fluctuations in the (small) 
number of pulsars in this study. 

6. Sources in the NVSS 

The NRAO VLA Sky Survey (NVSS, Condon et al. 1998) 
is a survey of the sky above declination —40° at 1400 MHz. 
Its final resolution is 45" and is comparible with the 
WENSS resolution of 54" (in right ascension). The band- 
width was effectively 42 MHz. The noise in the final NVSS 
maps is about 0.45 mJy. These maps are constructed by 
taking the average of a number of snapshot observations. 
Each point on the sky is observed in about three snap- 
shots. 

Seventeen of the twenty-five pulsars with a WENSS 
counte rpart are also detected in the NVSS (see Kaplan 
et al. fll998j ), who did not include PSR B0138+59 and 



Han & Tian Jl999| )). Two pulsars (PSR B0353+52 and 
PSR J1518+4904) with a marginal WENSS c ounte rpart 
coincide with a NVSS source. Kaplan et al. ( 1998 ) also 
lists the WENSS sources which coincide with their NVSS 
sources. Besides the three pulsars already mentioned, 
PSR B0809+74 is also not in their list, though it matches 
their coincidence criterion. Eight pulsars with a WENSS 
counterpart are not detected in the NVSS. This can be 
caused by scintillation or these pulsars must have a steep 
spectrum (spectral index between about —1.6 and —2.0). 
Two pulsars (PSRs J0218+4232 and B2319+60) may have 
a very steep spectrum (spectral index less than —2.6 and 
—2.3, respectively). 

From the WENSS and NVSS flux density values the 
spectral index can be calculated. These can be compared 
with the spectral indices as determined from a large num- 
ber of dedicated pulsar observations, as by LYLG. The 
latter indices will suffer less from scintillation effects. The 
uncertainty due to scintillation is included in the errors 
given by LYLG. 



-1 -0.5 0.5 1 

PSR Spectral Index - WENSS-NVSS Spectral Index 

Fig. 8. Histogram of the differences of the spectral index 
determined from the WENSS and NVSS source flux den- 
sities and of the pulsar spectral index, as determined by a 
large number of dedicated pulsar observations (Lorimer et 
al. 1995). The dashed line is a Gaussian distribution with 
the same mean and standard deviation as the histogram. 

From Table | it is clear that there are some large dif- 
ferences between the spectral index as determined by the 
WENSS and NVSS flux densities and the long term aver- 
aged spectral index. The differences are displayed in Fig. |[ 
This histogram has a mean of —0.14 and a standard devi- 
ation of 0.48. There is a large spread and a WENSS-NVSS 
source flux density comparison is not a very precise way to 
determine a spectral index. This will limit the success of a 
pulsar candidate selection based on their WENSS-NVSS 
spectral index. 

7. Discussion 

The deviations between the flux density values of pulsars 
and the correlated WENSS and NVSS sources are mainly 
caused by scintillation effects. Small-scale inhomogeneities 
in the interstellar medium affect the travel path of radio 
waves and can amplify or weaken them. For reviews, see 
Rickett qi99CD or Narayan ( |1992| ). Walker ( |1998| ) summa- 
rizes the involved equations and dependencies. 

In the case of strong scintillation, the phase changes 
due to the scattering in a certain region are larger than the 
changes in phase due to normal geometry. The scintillation 
is called weak in the opposite case. The boundary between 
the two is dependent on the frequency of the radio waves 
and the distance to the pulsar. At 325 and 1400 MHz 
almost all pulsars are in the strong scintillation regime. 

Two types of strong scintillation exist: diffractive and 
refractive. They differ in their typical timescale r, fre- 
quency bandwidth Ai/ and size of the resulting flux den- 
sity variation. This strength of the scintillation is usually 
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Fig. 9. Grayscale WENSS map of PSR B0329+54. The 
plot is centered around the pulsar position. To improve 
the visibility of the spokes, the greyscale range has been 
adjusted (top panel). The grey oval in the lower left corner 
indicates the FWHM beam size. The spokes are caused 
by diffractive scintillation, the ring is caused by refractive 
scintillation. 



quantified by the modulation index, i.e. the rms fractional 
flux density variation. 

In the case of strong scintillation waves from multiple 
locations in the scattering region interfere constructively 
(or destructively). Both the typical timescale and the fre- 
quency bandwidth are small. The modulation index equals 
one. The timescale is dependent on the relative velocities 
of the pulsar, the interstellar medium and the Earth. For 
the WENSS observations of pulsars typical timescales are 
1 to 10 minutes and typical bandwidths are 10 Hz to 500 
kHz. This means that for almost all pulsars any variations 
due to diffractive scintillation are averaged out over the 5 
MHz WENSS bandwidth and when the 6 x 18 short ob- 
servations spread over 6 x 12 hours are combined. 

Only for PSR B0809+74 are the diffractive scintilla- 
tion timescale and bandwidth large enough that some ef- 



fect remains. From the equations given by Walker (1998) 



one finds a timescale of 12 minutes and a bandwidth 
of 500 kHz. The actual value of the typical scintillation 
timescale and bandwidth are even higher, since several au- 



thors have already shown that the Taylor-Cordes distance 
model ( [Taylor fc Cordes, 1993| ) gives too small predictions 
for this pulsar (e.g. Rickett et al. 2000). Diffractive scintil- 
lation can explain the WENSS image of this pulsar (Fig. 
||) . As 40 minutes is the time between two observations of 
a field in one mosaic, flux density variations on that time 
scale cause spokes in the map. Such spokes can also be 
seen in the complete map of PSR B0329+54 (Fig. |). 

Refractive scintillation is caused by the focussing ef- 
fect of a large scattering region. The timescales and band- 
widths involved are much larger than in the case of diffrac- 
tive scintillation. The modulation index is also smaller. 
The scintillation bandwidth is of the order of the observing 
frequency. Refractive time scales for the pulsars detected 
in the WENSS vary from a couple of days to several years 
and the expected modulation index from 0.05 to 0.3. 

If the refractive time scale is less than the time be- 
tween observations of the same mosaic (couple of days to 
several years), any flux density variation will be averaged 
out when the mosaics are combined. However, stong flux 
density variations between 12 hour sessions cause a ring 
at the first grating ring of the synthesised beam. Since the 
mutual distances between the dishes are multiples of 72 m, 
the ring will have radius of 72 m / (c/325 MHz) radians, 
i.e. 44' in right ascension and 44' x cosec S in declination. 
The second and higher grating rings are not visible, since 
data far from the field center gets a low weighting factor 
when the final image is created. 

It is hard to give good estimates for the expected 
diffractive and refractive scintillation timescales. They de- 
pend on the often poorly known pulsar velocity. For large 
pulsar velocities (compared to the velocity of the interstel- 
lar medium, being about 50 km/s) the dependence is as 
one over the square root of this velocity. Since some pul- 
sar velocities might be up to several hundred kilometers 
per second, this cannot be neglected. I divided the pulsars 
that are detected in the WENSS in two groups, based on 
their expected refractive modulation index if their velocity 
is neglected. Both groups had similar relative deviations 
between their measured and expected flux densities. 

The observed modulation index is about 0.4, much 
larger than the expected value of about 0.2. The WENSS 
pulsar flux densities vary more due to refractive scintilla- 
tion than predicted by the equations. This has been ob- 
served before by several authors and is attributed to the 
assumption that the turbulence in the interstellar medium 
has a Kolmogorov spectrum (e.g. Blandford & Narayan, 
1985). 

The NVSS flux densities are even more affected by 
scintillation effects. Each point in the NVSS maps is an 
average of about three snapshots. Two of these three are 
taken right after each other (snapshot series are taken at 
constant declination and increasing right ascension, see 
figure 7 in Condon et al. 1998) and very little averaging 
takes place. The expected refractive modulation index at 
1400 MHz is larger than at the WENSS frequency. The ex- 
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pected diffractive frequency bandwidth is also larger at the 
NVSS frequency, even relative to the increased total band- 
width used in the NVSS. One therefore expects the differ- 
ences between the pulsar flux densities in the NVSS and 
the flux densities reported by LYLG to be larger than the 
differences reported in this study. This is indeed observed: 
the modulation index of the WENSS sources in Table |, 
excluding PSRs J0218+4232, B0329+54 and B2021+5f is 
0.40, the modulat ion in dex of the NVSS sources in table 
1 of Han & Tian (|l999|) is .60. 

De Breuck et al. ( 2000|) show that the total spectral 
index distribution of WENSS sources that are correlated 
with a NVSS source (5 I 32 5 > 50 mJy and \b\> 15°) has 
a mean of —0.80 and a standard deviation of 0.24. Pul- 
sar have a mean spectral index of —1.6 (see LYLG). An 
increase of 0.5 due to scintillation will move the pulsar 
spectral index well into the distribution of normal sources, 
which are much more frequent. The spectrum of some 
point-like quasars will also be affected by scintillation and 
some of them will seem to have a much steeper spectrum 
than they really have. This effect should be taken into con- 
sideration if pulsar candidates are selected on the basis of 
their spectral index derived from the WENSS and NVSS. 
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Table 4. List of the pulsars with no WENSS counter- 
part, their expected flux densities at 325 MHz and the 
local noise in the WENSS map. The expected flux den- 
sity for pulsars B0841+80, J1012+5307 and B1639+36A 
is equal to their 400 MHz flux density. The flux density 
of PSR J2044+46 is not known. 'ND' in the last column 
indicates, that the expected pulsar flux density is higher 
than three times the local noise level. 



DQD 

name 


PSR 5325 

estimate 
(mJy) 


Q 

>J noise 

(mJy) 




nnm i i A T 


19 ± 


5 


Q O 
O.O 


"NT "Hi 


JdUU4oH-oo 


> 3.5 ± 


0.5 


O.O 






4.1 ± 


0.5 


o n 

o.y 




TiAncn i fie 


11.4 ± 


1.3 


6.1 


1NJJ 


B0105+65 


14.6 ± 


1.9 


3.2 


ND 


T2H1 A A i CO 


7.0 ± 


0.5 


Q A 

o.4 




"Dnoot; I 7H 
dUzzo-\- (U 


3.5 ± 


0.4 


Q C 
O.O 




"DriQQ i i a c 
r5Uoo±+4o 


14 ± 


3 


A O 

4.z 


"NT "Pi 

IN JJ 


DUooy+Oo 


4.9 ± 


1.0 


a n 
4.U 




B0410+69 


10.5 ± 


0.7 


4.2 




Drift a o i on 
r3U643+oO 


11.2 ± 


1.8 
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"NT "T"\ 

INJJ 


rJU/Ol+OZ 


12.4 ± 


2.5 


O.O 


INJJ 


rJUo41+oU 


1.5 




O 1 




rJUyl / +63 


6.9 ± 


1.4 


Q C 
O.O 




J1012+5307 


30 ± 30 


4.2 


ND 


T41 OOP 1 Q9 


15 ± 


7 


O 

z.z 


IN D 


-DlOoa+oDA 


3.0 




A Q 

4.o 




rJl /o3+o2 


5.1 ± 


1.1 


3.5 




B1811+40 


18.4 ± 


1.0 


3.3 


ND 


B2000+32 


8.3 ± 


0.5 


7.3 




B2002+31 


27.9 ± 


1.8 


7.8 


ND 


B2011+38 


34.3 ± 


1.8 


8.5 


ND 


B2022+50 


7.1 ± 


0.8 


4.8 




B2027+37 


27 ± 


3 


7.4 


ND 


B2035+36 


8.7 ± 


2.2 


7.0 




B2036+53 


4.6 ± 


0.8 


5.4 




J2044+46 






6.4 


? 


B2045+56 


6.6 ± 


1.8 


4.8 




B2241+69 


3.6 ± 


1.4 


4.2 




B2303+30 


38 ± 


4 


3.7 


ND 


B2303+46 


1.7 ± 


0.7 


3.8 




B2334+61 


16 ± 


4 


6.9 




B2351+61 


23 ± 


4 


5.5 


ND 



Table 5. List of spectral indices from Lorimer et al. (1995) 
and from a WENSS-NVSS flux density comparison for all 
pulsars with both a WENSS and a NVSS counterpart. 
The spectral index for PSR J1518+4904 was taken from 
Kramer et al. (1999). The errors on the pulsar spectral 
index (second column) include variations due to scintilla- 
tion, the errors on the WENSS-NVSS spectral index do 
not include scintillation effects. 
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Spec.ind. 


Spec.ind. 


name 
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